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ON THE ORIGIN OF THE COLOR-MAGNITUDE RELATION IN THE VIRGO CLUSTER 
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ABSTRACT 

We explore the origin of the color-magnitude relation (CMR) of early type galaxies in the Virgo cluster 
using spectra of very high S/N ratio for six elliptical galaxies selected along the CMR. The data are 
analysed using a new evolutionary stellar population synthesis model to generate galaxy spectra at the 
resolution given by their velocity dispersions. In particular we use a new age indicator that is virtually 
free of the effects of metallicity. We find that the luminosity weighted mean ages of Virgo ellipticals 
are greater than 8 Gyr, and show no clear trend with galaxy luminosity. We also find a positive 
correlation of metallicity with luminosity, color and velocity dispersion. We conclude that the CMR 
is driven primarily by a luminosity-metallicity correlation. However, not all elements increase equally 
with the total metallicity and we speculate that the CMR may be driven by both a total metallicity 
increase and by a systematic departure from solar abundance ratios of some elements along the CMR. 
A full understanding of the role played by the total metallicity, abundance ratios and age in generating 
the CMR requires the analysis of spectra of very high quality, such as those reported here, for a larger 
number of galaxies in Virgo and other clusters. 

Subject headings: galaxies: abundances — galaxies: clusters: general — galaxies: clusters: individual 
(Virgo) — galaxies: elliptical and lenticular, cD — galaxies: evolution — galaxies: 
stellar content 



1. INTRODUCTION 

To understand how galaxies form and evolve, we need to 
account for the color- magnitude relation (CMR): luminous 
early-type galaxies in clusters are observed to be redder 
than fainter ones (Visvanathan & Sandage 1977; Bower et 
al. 1992a). This tight relation, together with the related 
Mg2-cr relation (CoUess et al. 1999) links the mass of a 
galaxy, through its luminosity, to its constituent stellar 
populations. 

The origin of the CMR is hotly debated: it could be 
caused either by a variation of the mean stellar metallicity 
along the sequence (Arimoto & Yoshii 1987; Kodama & 
Arimoto 1997), or by a combination of age and metallicity 
variations (Gonzalez 1993; Ferreras, Chariot & Silk 1999). 
Larson (1974) suggested a dissipational collapse picture for 
the formation of ellipticals with an early formation of the 
bulk of the stars. In this scenario more luminous galax- 
ies have greater binding energies and can therefore achieve 
higher metallicitics. Tcrlevich et al. (1999) provided sup- 
port for this view by showing that the Coma cluster galax- 
ies on the CMR are old, and that the galaxies with bluer 
colors than expected for their luminosity are younger. In 
more complex star formation histories, such as implied by 
hierarchical merging (Cole et al. 1994; Kauffmann & Char- 
lot 1998) where galaxies have some recent star formation 
due to late merger events, the CMR might be expected to 
arise from a combination of age and metallicity variations. 
Indeed Gonzalez (1993), J0rgensen (1999) and Trager et al. 



(2000a) find strong evidence for a significant intermediate 
age population in some elliptical galaxies. Trager et al. 
(2000b) reported a large age range in their full sample but 
note that the galaxies in clusters are generally old. They 
also find a tight age-metallcity-velocity dispersion relation 
such that at a fixed velocity dispersion, metal rich galaxies 
are young and conjecture that this relation might give rise 
to a tight CMR. 

Broad band colors and line strength indices of integrated 
stellar populations are sensitive to both changes in age 
and metallicity and this degeneracy has prevented an un- 
ambiguous analysis of the ages and metallicitics of ellip- 
ticals along the CMR (Worthey 1994). Even the widely 
used age indicator in the Lick/IDS system, H/?, shows a 
non-negligible dependence on metallicity (Worthey 1994) 
and can be affected by emission (Davies, Sadler & Peletier 
1993; Gonzalez 1993). The origin of the CMR is not yet 
fully understood, in part, because of the lack of appropri- 
ate tools of analysis. 

In this Letter we study the stellar populations of six 
Virgo ellipticals selected along the CMR for which we have 
very high S/N spectra. We use a new stellar population 
synthesis model in combination with new, more sensitive 
age indicators to derive age estimates that are virtually 
independent of the metallicity. 

2. STELLAR POPULATION MODELS 

Here we use the evolutionary stellar population syn- 
thesis model of Vazdekis (1999), which employs the ex- 
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tensive empirical stellar spectral library of Jones (1999). 
The model predicts spectral energy distributions (SEDs) in 
the optical wavelength range for simple stellar populations 
(SSPs) at a resolution of 1.8 A (FWHM). Previous mod- 
els (e.g., Worthcy 1994; Vazdekis et al. 1996) used mostly 
the Lick polynomial fitting functions (Worthey et al. 1994; 
Worthey & Ottaviani 1997) to relate the strengths of se- 
lected absorption features to stellar atmospheric param- 
eters. The fitting functions are based on the Lick/IDS 
stellar library (FWHM~9 A, Worthey et al. 1994), thus 
limiting the application of the models to strong features. 
Furthermore the Lick stellar library has not been flux cal- 
ibrated, so offsets had to be applied to data obtained with 
different spectrographs (Worthey & Ottaviani 1997). 

The new model provides flux calibrated, high resolution 
spectra, therefore allowing an analysis of galaxy spectra 
at the natural resolution given by their internal velocity 
broadening and instrumental resolution. Furthermore, as 
the model outputs are SEDs rather than predicted index 
strengths, it is easy to define new indices and even to con- 
front the model predictions with the detailed structure of 
observed absorption features. Vazdekis & Arimoto (1999) 
defined a new age indicator [H7-I- i(FeI-fMgi)]cr (hereafter 
H7o-), which is virtually free from metallicity dependence. 
To cover a large range in cr, Vazdekis & Arimoto defined 
three slightly different indices which give stable and sen- 
sitive age predictions within the a ranges quoted in the 
sub-indices: H7ioo<ct<i75, H7i5o<cr<225 and H7225<^<300- 

We emphasise that the model spectra we are fitting have 
a single metallicity and age. We are measuring luminosity- 
weighted properties, which means that the presence of a 
young (i.e., bright) stellar component will disproportion- 
ally affect the mean age and therefore will be detected 
more easily. 

3. DATA 

We obtained major axis spectra of very high S/N ratios 
for six Virgo ellipticals NGC 4365, NGC 4387, NGC 4464, 
NGC 4473, NGC 4478, and NGC 4621, selected along 
the CMR of Bower, Lucey & Ellis (1992a,b) (see Fig. 
and also Table 1). The observations were performed at 
the WHT (4.2m) at the Observatorio del Roque de Los 
Muchachos, La Palma, in April 21-22, 1999. We used the 




Fig. 1. — The CMR in the Virgo cluster. Data and a best-fit 
regression (ellipticals only) from Bower, Lucey & Ellis (1992a, b). 
For NGC 4464 we have taken Vt from RC3. For NGC 4478 we 
used Michard (1982) U — V color (60" aperture) applying an offset of 
—0.047, obtained by a comparison of his non dereddened U—V colors 
with those of Bower, Lucey & Ellis (1992b) which are dereddened by 
—0.02. Solid symbols represent our sample of galaxies, the shapes 
are used to identify the six galaxies in the following figures. 



ISIS double-beam spectrograph with the EEV12 CCD and 
R600B grating in the blue channel, which provided a dis- 
persion of 0.44 A/pix and a spectral range of AA ~4000- 
5500 A. The slit width was set to 1'.'6 giving a spectral 
resolution of 2.4 A (FWHM). The seeing was 2-3" in the 
first night and ~1'.'5 in the second. The CCD was binned 
x3 in the spatial direction giving 0'.'57 per bin. Multiple 
exposures of 35 min length were obtained for each galaxy. 

Data reduction was done with standard IRAF packages. 
We made use of tungsten flatfields obtained for each galaxy 
pointing. The continuum shape was corrected to a relative 
flux scale with spectrophotometric standards. 

We summed up the spectra within rg/lG, where effective 
radii were taken from Burstein et al. (1987). Before sum- 
mation the central galaxy spectra were carefully aligned 
with the rows of the array, and each row in spatial direc- 
tion was corrected for the velocity shift due to rotation. 
Variations of ct within re/ 10 were small (< 25 km s~^) for 
all galaxies and therefore no correction was applied. Ta- 
ble 1 hsts the S/N ratio per A achieved in the summed 
spectra around the H7 feature. 

4. THE AGES AND METALLICITIES OF VIRGO 
ELLIPTICALS 

In this section (i) we compare the sensitivity and preci- 
sion of different age indicators, (ii) demonstrate that not 
all metals are enhanced at the same rate in more luminous 
(higher dispersion) galaxies and (iii) show that for the six 
Virgo ellipticals studied here metallicity increases with lu- 
minosity, velocity dispersion and color but no trend with 
age is apparent. 

Our primary goal is to measure the ages of the central 
regions of the galaxy sample along the CMR. We use dif- 
ferent age indicators: Hja, H^f and II/3 (as defined by 
Vazdekis & Arimoto (1999), Worthey & Ottaviani 1997, 
and Worthey et al. 1994, respectively). Instead of ap- 
plying a corrections to the line strength measurements we 
modified the model predictions. We divided our sample 
into three groups containing galaxies with similar velocity 
dispersions and smoothed the SEDs of Vazdekis (1999) ac- 
cordingly. We chose (a) atot = 145, (b) atot = 190 and (c) 
atot = 260 kms-i (where cr^^t = cr^aia^y + af^^^^ ) . To min- 
imize the effect of small velocity dispersion differences the 
spectrum of NGC 4387 was convolved with a Gaussian of 
cr = 85 kms-i and that of NGC 4621 with cr = 105 kms~^ 

Fig. I shows age/metallicity diagnostic diagrams for the 
age indicators as a function of the mean metallicity indica- 
tor [MgFe] (Gonzalez 1993). Velocity dispersion increases 
in each column from top to bottom. Note that for a given 
age and metallicity the strength of the [MgFe] index is 
smaller for larger a. The left panels of Fig. || show how 
well H7o- is able to break the age-metallicity degeneracy, 
much better than, e.g., HSp. The lines of constant age 
are essentially horizontal, which means that a given mea- 
surement of Il7cr corresponds to a unique age. H/3 is also 
a good age indicator, the new models predict that H/3 is 
more sensitive to age, and less sensitive to metallicity than 
models using the Lick fitting functions. 

Overall we find good agreement between the ages in- 
ferred from H7o- and II/3. All our galaxies, except for 
NGC 4478, are older than - 11 Gyr. NGC 4478 is 8-9 Gyr 
old and has aU — V color 0.06 mag bluer than other Virgo 
ellipticals of equivalent luminosities (see Fig. 
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Fig. 2. — Plots of H7o-, U5p and H/3 (from left to right) versus the 
metallicity index [MgFe] (defined in Gonzalez 1993). The central ve- 
locity dispersion increases in three groups from top to bottom, rep- 
resenting the CMR. Overplotted are the models by Vazdekis (1999). 
Lines of constant [Fe/H] = —0.7, —0.4, 0.0 and -fO.2 are shown by 
thick dot-dashed, thick dotted, thick solid and thick dashed lines, 
respectively. Thin dotted lines represent models of constant ages 
which are quoted in Gyr. 

The middle panels of Fig. ^ show that the age estimates 
based on HSp remain strongly coupled with metallicity 
estimates. (NB. the Hjf index is similar and is not plot- 
ted.) The age estimates based on RSp are younger than 
those derived from Hj^ or H/3 but we are unable to dis- 
tinguish between genuinely younger ages and the possible 
systematic effects of non-solar abundance ratios. For ex- 
ample, using HSp or H7F as the age indicator one derives 
younger ages when Mgb and/or Mg2 are used as the metal 
indicator than when the Fe lines alone are used (Worthey 
1998; Kuntschner 2000). 

Fig. 1^ plots the H'ja indices against the strength of four 
metal absorption lines and illustrates the effects of non- 
solar abundance ratios in determining metallicity. The 
left-hand columns in Fig. ^ show only a very modest in- 
crease in the metallicity derived from Fe3|^ or Ca4227 as 
velocity dispersion increases from top to bottom. We note 
that the G-band also follows Fe3. In contrast the right- 
hand columns in Fig. |^ show that the metallicity derived 
from Mg b or CN2 increases significantly towards high ve- 
locity dispersion galaxies, extending beyond the model 
grids in some cases. Although Ca, like Mg, is an a-element, 
we find that Ca4227 does not track Mgb confirming the 
results reported by Vazdekis et al. (1997), Worthey (1998) 
and Trager et al. (2000a). We conclude that there are at 
least two families of metal lines which give rise to very 
different metallicity estimates and exhibit different trends 
along the CMR. Despite these complexities in estimating 

lFe3 = i(Fe4383 -I- Fe5270 + Fe5335), Kuntschner (2000). 




Fe3 Ca4227 Mgb CN2 

Fig. 3.— Plots of vs Fe3, MgbjjCa4227 and CN2. Symbols 
and line styles are the same as in Fig. eI 

"metallicity" , we emphasize that by using Hja these un- 
certainies do not influence our derived ages. 

Fig. ^ shows our estimates for luminosity weighted mean 
age and metallicity (derived from 117^ and H/3 vs. [MgFe] 
in Fig. 1^) versus galaxy luminosity, color, and velocity dis- 
persion. Except for NGC 4478, all ellipticals are old with 
ages ranging from ~11 Gyr to ~20 Gyr. There is no trend 
in age along the CMR but there is a clear, positive corre- 
lation of metallicity with V^-, U — V and a. If we correct 
the metallicity estimates to the value they would have if 
all the galaxies were 14 Gyr old the scatter in Fig. ^ is 
reduced. Furthermore our [MgFe] index yields very simi- 
lar metallicities to those derived by correcting the Fe in- 
dices for non-solar abundance ratios, following Trager et 
al. (2000a). We conclude that for these six galaxies in 
Virgo the metal content increases along the GMR. 

5. CONCLUSIONS 

We have performed a detailed spectral analysis of six 
elliptical galaxies in the Virgo cluster. We have combined 
very high S/N spectra with a new evolutionary stellar pop- 
ulation synthesis model. This model enables us to study 
galaxy spectra at the resolution given by their velocity dis- 
persions and determine ages that are virtually free of any 
dependence on metallicity. We find that elliptical galaxies 
in Virgo are generally old. 

We find that the metallicity as measured by [MgFe] cor- 
relates with luminosity, U — V and velocity dispersion. 
We conclude that the CMR in Virgo is driven primarily 
by an increase in metallicity with luminosity. We find that 
not all elements increase equally with the total metallic- 
ity estimates. For example the increase in metallicity as 
measured from the Fe or Ca indices is very modest along 
the CMR, whereas the Mg or CN indices indicate a more 
rapid increase with galaxy luminosity. This is consistent 
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Fig. 4. — Mean galaxy ages and metallicities as estimated from 
H7ct (filled symbols) and H/3 (open symbols) (derived from Hfl vs 
[MgFe] diagrams) vs V!^ , U — V and velocity dispersion. The error 
bars for the H7ct age and metallicity estimates are mainly due to 
noise, while those for H/3 are dominated by the larger age/metallicity 
degeneracy affecting this index. 



with the findings of Worthcy, Fabcr, & Gonzalez (1992), 
Kuntschncr (2000) and Tragcr ct al. (2000b) who find an 
increase of the [Mg/Fc] ratio with increasing central veloc- 
ity dispersion. We speculate that the CMR in Virgo may 
be driven by both a total metallicity increase and by a 
systematic departure from solar abundance ratios of some 
elements along the CMR. 

The luminosity weighted mean ages range from ~ 8 — 
20 Gyr but there is no clear trend between age and V^^, 
U - V or a. Trager et al. (2000b) note that at a 
given velocity dispersion galaxies with younger ages have 
higher metallicity. Our Virgo galaxies fall on the Trager 
dispersion-age-metallicity relation but span less than half 
the range of ages and metallcities that are found in their 
sample. Trager ct al. suggest that the distribution of ages 
and metallicities amongst ellipticals may depend on en- 
vironment and indeed our Virgo galaxies exhibit a much 
smaller range than that found in their mixed environment 
sample. We cannot determine whether the small age range 
wc find here is universal in cluster ellipticals but we intend 
to address this question in future work. 
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Table 1 
Virgo galaxy sample. 



NGC 


Exposure (minutes) 


re/10 


S/N per A 


a{kms ^) 


V° ^ 


{U-V)% b 


4464 


105 


0'.'6 


~160 


~135 


12.46"= 


1.46* 


4387 


105 


2'.'6 


~130 


~105 


11.99 


1.50 


4478 


105 


1'.'7 


~135 


~135 


11.23 


1.50 


4473 


70 


2'.'7 


~280 


~180 


10.06 


1.61 


4621 


140 


6'.'0 


~490 


~230 


9.81 


1.64 


4365 


70 


6'.'7 


~250 


~255 


9.62 


1.67 



''Bower, Luccy & Ellis (1992b). 

'^Michard (1982). {U — corrected to re/10 from {U — V)e using the mean 
gradient of Trager et al. (2000b). Galactic extinction corrected using Schlegel et 
al. (1998). 

^RC3. 



